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Abstract 


We set up a formalism for calculating transverse-momentum-dependent parton distribution func¬ 
tions (TMDs) of a large nucleus using the tools of saturation physics. By generalizing the quasi- 
classical Glauber-Gribov-Mueller/McLerran-Venugopalan approximation to allow for the possibil¬ 
ity of spin-orbit coupling, we show how any TMD can be calculated in the saturation framework. 
This can also be applied to the TMDs of a proton by modeling it as a large “nucleus.” To illustrate 
our technique, we calculate the quark TMDs of an unpolarized nucleus at large-x: the unpolarized 
quark distribution and the quark Boer-Mulders distribution. We observe that spin-orbit coupling 
leads to mixing between different TMDs of the nucleus and of the nucleons. We then consider 
the evolution of TMDs: at large-x, in the double-logarithmic approximation, we obtain the Su- 
dakov form factor. At small-x the evolution of unpolarized-target quark TMDs is governed by 
BK/JIMWLK evolution, while the small-x evolution of polarized-target quark TMDs appears to 
be dominated by the QCD Reggeon. 
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I. INTRODUCTION 


Over the past decade quark and gluon transverse momentum-dependent parton distribu¬ 
tion functions (TMDs) jl, 2] have become an integral component of our understanding of 
the momentum-space structure of the nucleon. At the same time the main principles for 
calculating TMDs in the perturbative QCD framework have remained essentially the same 
over the years: one parameterizes the initial conditions at some initial virtuality Q 2 = Ql 
and then applies the Collins-Soper-Sterman (CSS) evolution equation lj] to find the TMDs 
at all Q 2 . The initial conditions are non-perturbative, and have to be constructed using 
models of the non-perturbative QCD dynamics (see 3] and references therein). Since very 
little is known about non-perturbative effects in QCD, often one uses the same form for 
the parameterization of the initial conditions for several different TMDs, frequently assum¬ 
ing a Gaussian dependence of the TMDs on the parton transverse momentum [ 4 ]. It is 
clearly desirable to have a better control of our qualitative and quantitative understanding 
of TMDs. 

To this end one can employ the recent progress in our understanding of small-a; physics 
and parton saturation [5| [l_l| to put constraints on the TMDs and even calculate them in 
the high-energy limit. When calculating TMDs one often works either in the s rs./ Q 2 » k 2 
(large-x) or in the s Q 2 k^ (small-x) regimes, where s is the center-of-mass energy and 
x = Q 2 /(s + Q 2 ) if one neglects the proton mass. I11 either case the energy s is large, and the 
techniques of high-energy QCD should apply. The degrees of freedom in saturation physics 
are infinite Wilson lines along (almost) light-like paths. The definition of the TMDs involves 
light-cone Wilson lines as well jl2, 13], although the integration paths are semi-infinite, 
forming the so-called “light-cone staple”. We can see that there are both similarities and 
differences between saturation physics and the physics of TMDs. The interface of these two 
sub-fields of quantum chromodynamics (QCD) has been explored in [14l-|24]. 

In the past some_success has been achieved in applying saturation physics to study the 
Sivers function 
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26| both in semi-inclusive deep inelastic scattering (SIDIS) and in the 


Droll-Yan process (DY). In |27J the Sivers function was constructed by generalizing the 


quasi-classical Glauber-Gribov-Mueller (GGM) [28] /McLerran-Venugopalan (MV) 29 


3ll] approximation of a heavy nucleus with atomic number A 1. The presence of the 
atomic number generates a resummation parameter o? s A 1 / 3 [32, |33j allowing a systematic 
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resummation of multiple rescatterings, which are essential for the Sivers function. This 
picture can also be applied to the proton if one models it as a large “nucleus.” This large- 
nucleus approximation is known to work well in describing the data from deep inelastic 
scattering (DIS) experiments on a proton at low-a;; it is therefore possible that it would give 
a reasonable description for proton TMDs as well. The result of 27j was an explicit form of 
the Sivers function in the s ~ Q 2 regime, which was different from a simple Gaussian 

in kx and which can be used as the initial condition for its Collins-Soper-Sterman (CSS) 
evolution [1], Another important result of was the realization that the Sivers function 
can be produced via two different channels: one is the standard “lensing” mechanism of 
with additional momentum broadening due to multiple rescatterings in the nucleus, 
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while the other mechanism was due to the orbital angular momentum (OAM) of the nucleus 
combined with multiple rescatterings. This latter channel has not been reported before [2^], 
and it dominates mainly in the regime where multiple rescatterings are important, or, more 
precisely, for kx < Q s /\/<Ts, where Q s is the saturation scale of the nucleus. It appears that 
applications of saturation physics to the calculation of TMDs may lead to qualitatively new 
channels of generating the relevant observables. 


The aim of this work is twofold. First of all we want to generalize the approach of [21 
to the calculation of any TMD in the quasi-classical approximation and for s ~ Q 2 kx- 
This is accomplished in Sec. [TT| for the case of an unpolarized nucleus; the generalization to 
the polarized case is straightforward and is left for future work. The quasi-classical TMD 
calculation is accomplished using the factorization given in Eq. (TT8|) (and, in more detail in 
Eq. (l39li ). which is constructed by analogy with the particle production cross section: the 
TMD is a convolution of the classical Wigner distribution W describing the nucleons in the 
nucleus with the TMD of one of the nucleons (4> N ) and the Wilson lines S of the “staple” 
evaluated in the GGM/MV approximation. The functional form of the nucleon Wigner 
distribution in the nucleus can be constrained using the symmetries of the nuclear ground 
state. Considering for simplicity an unpolarized nucleus we arrive at the parameterization of 
Eq. (J55j). With the help of this unpolarized nuclear Wigner distribution we construct the un¬ 
polarized quark distribution in Eq. (16TT) and the Boer-Mulders distribution 37] in Eq. (l78lh 
Note that the nuclear effects on the TMDs are not limited to the fc^-broadening as is com¬ 


monly assumed in the literature 


Other TMDs for the (transversely or longitudinally) 
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polarized nucleus and for the unpolarized nucleus could be constructed by analogy, but we 
limit the examples to the two TMDs mentioned above along with the Sivers distribution 
constructed in [27 ]. 

An interesting consequence of the calculation is a mixing between different TMDs of the 
nucleus and of the nucleons. For instance, as was already observed in \l7}, the nuclear Sivers 
function is a linear superposition of the nucleon Sivers function and the unpolarized quark 
distribution of the nucleon. Similarly, we find that the unpolarized quark distribution of a 
nucleus in Eq. (IbTf is a linear superposition of contributions due to the unpolarized nucleon 
quark distribution and the Sivers function of the nucleons. The nuclear Boer-Mulders func¬ 
tion in Eq. (I75|) receives contributions from the nucleon Boer-Mulders function, transversity 
and pretzelosity. Such TMD mixing means that, in the quasi-classical approximation, dif¬ 
ferent nuclear TMDs are expressed in terms of the same Wigner distributions and various 
nucleon TMDs. This may lead to relations between different nuclear TMDs which, in prin¬ 
ciple, could be tested experimentally. 


The second aim of this work is to understand TMD evolution using the methods of 
saturation physics. In Sec. IIIII we consider TMD evolution at large and small values of 
Bjorken-x. The large-x regime, corresponding to s ~ Q 2 k^, is considered in Sec. IIII.ll 
There we show that the TMD evolution can be thought of as the evolution of the light-cone 
Wilson lines in the “staple” (if one works in the light-cone gauge of the probe), as shown in 
Fig. [6] In the large-x regime the Wilson lines combine to form an operator corresponding to 
the propagation of a (quark or gluon) dipole from the nucleus to positive light-cone infinity 
(S1DIS) or from negative infinity to the nucleus (DY). (A fundamental semi-infinite dipole is 
used for quark TMDs, while a semi-infinite adjoint gluon dipole should be used for the gluon 
TMDs.) In the double-logarithmic approximation resumming powers of a s In 2 (Q 2 /k^) we 


recover the standard Sudakov form factor 


which one also obtains from the CSS evolution. 


In the low-x regime (s ^ Q 2 fcy) the diagrams giving the leading contribution to 
the TMDs are slightly different from the large-x graphs, as illustrated e.g. in Fig. [HI For 
unpolarized-target quark TMDs in the leading single-logarithmic approximation, which re¬ 
sums powers of a s ln(l/x), we show in Sec. 1111.21 that TMD evolution is driven by the 
evolution of a fundamental dipole scattering amplitude on a nucleus, that is by the evo¬ 
lution of the dipole made out of infinite Wilson lines. The evolution of this object is 
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we. 
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1-known: in the large- N c limit it is given by the Balitsky-Kovchegov evolution equation 


.431 ], while beyond that limit it can be found using the Jalilian-Marian-Iancu-McLerran- 


Weigert-Leonidov-Kovner (JIMWLK) evolution equation 44j-|47j. For the polarized-target 
TMDs the situation is more subtle, as the evolution has to carry the polarization infor¬ 
mation. We argue that polarized quark TMDs at small -x are governed by the nonlinear 


version 


481 ] of the QCD Reggeon evolution equations 49-1551] (see Fig. [T3l and Eq. (j94h). 


which employs the dipole amplitude obtained from the BK/JIMWLK equations. 

We have thus constructed quasi-classical expressions for TMDs at large-.!' in the GGM/MV 
approximation, and constructed evolution of TMDs at both large- and small-.!. The results 
of this work can be used for constructing a global TMD fit. 


II. QUASI-CLASSICAL INITIAL CONDITIONS 


In this Section we expand and generalize the approach developed in [27]. While the 


quasi-classical method detailed below was used in 


27] to calculate the Sivers function of a 


large nucleus, here we will apply it to calculating TMDs of an unpolarized nucleus, such as 
the unpolarized quark distribution and the Boer-Mulders function. 


II. 1 . Quasi-Classical Factorization of the TMD Quark Distribution of a Heavy 
Nucleus 

II. 1.1. The TMD Decomposition 

The transverse-momentum-dependent (TMD) quark correlation function in a hadronic 
state | h(P, S)) with momentum P and spin S is defined by 

P,S) = |^=; J d 2 ~r e ik ' r (h(P,S)\ ?„(0)M[0.r]^(r) \h(P,S)) rt _„ (1) 

where a,/3 are Dirac indices, the separation vector between the quark fields is = 
(0 + ,r _ ,r), and the “staple-shaped” gauge link 7/[0,r] extends to future/past light-cone 
infinity (±oo - ) depending on the process under consideration. The abbreviated notation 
for the integration measure is d 2 ~r = d 2 r±dr~. We work in light-front coordinates 

!± = (x° ± ! 3 ) (2) 
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where g + ~ = 1 and g + ~ = 2 are two common choices of the light-front metric. We will 
use a frame in which the nucleus or the proton have a momentum predominantly in the 
x + -direction. 

Projecting the TMD correlator ([TJ) onto various Dirac matrices gives the distribution 
of quarks with zero, longitudinal, or transverse polarizations, and the independent spin- 
spin and spin-orbit correlations are parameterized in terms of boost-invariant TMD parton 
distribution functions (the “TMDs”). At leading power in the large momentum P + , the 


TMD decomposition can be written as (see, e.g. [371. l56l. l57j|) 
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where (S_, Sl) denote the transverse and longitudinal components of the spin vector, m is the 
hadron mass, and kxS_ = k x S y — k y S x with (x, y) the coordinates in the plane transverse 
to the beam. The eight leading-twist quark TMDs shown here 


{/i, fir , 9i, .9 it , h\ T , hi L , hi T , hi} 


( 4 ) 


are functions of x and kx = \k\. 

At lowest order in a s such that the gauge link W[0, r] can be neglected, the correlator (JT]) 
has a simple interpretation as the expectation value of the quark density in the hadronic 
state | h(P,S)): 


Mx,k, p, s) L = ^ (h(p, s)\bl,b k „ | h(p,s)) [(^(4))^ (u,(k)) 


(5) 


where D = g .j_ p + f d 2 b is a boost-invariant volume factor which normalizes the plane-wave 

states | h(P,S)). The gauge link W[0,r] is needed to make the definition gauge invariant, 


and it reflects t 
gauge fields 


re distortion of the quark distribution in the presence of initial- or final-state 


59 |] 


By tracing (J5]) with one of T G {y -1- , 7 + 7 5 ,7 5 7 + 7 j_}, we P r °j ec t out the distributions cor¬ 
responding to unpolarized quarks (U), longitudinally-polarized quarks (L), and transversely- 
polarized quarks in direction j (labeled T J ), respectively. This can be explicitly verified by 
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evaluating the associated spinor products in ([5]) using the spinor conventions of [601 ]: 


U: U a/ (k)-f + U a {k) = 2k + 

L : U a ,(k) 1 + 1 b U a (k) = 2k + 

T j : U a ,(k) 1 b 1 + 1 {U a (k) = 2k + 


1 0 
0 1 


cr'cr 


2 k + ill 


1 

0 


0 

-1 


G* <J 


2* + h 3 L 


0 5’ 1 - ii’ 2 


S^ 1 + i& 2 0 


cr'cr 


= 2 k + (S jl [a 1 ] 


+ <p ' 2 h 2 L„)- 


( 6 ) 

( 7 ) 

( 8 ) 


These projections select out the corresponding set of TMDs from the decomposition 


U : 0h + ] = iiv [07+] =fi~ 

2 L J m 

L : 0 [7+7S] = [0 7 + 7 5 ] = S L g 1 + ^ g 1T 

2 L J m 


( 9 ) 


T j : = ^Tr [0 7 5 7 + 7l] = S{h 1T + ^ S L hi L + ^^-=lhj; T + 4^, 


k{ 

m 


k{ (k ■ S) 


m m 


Li 

: ft, j 
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so that f\ represents the azimuthally symmetric distribution of unpolarized quarks in an un¬ 
polarized hadron, the Sivers function represents the azimuthally asymmetric distribution 
of unpolarized quarks in a transversely polarized hadron, and so on. 


II. 1.2. Quasi-Classical Factorization in a Heavy Nucleus 


As derived in 


271] . a heavy nucleus with A 1 nucleons in the quasi-classical approxi¬ 


mation a 2 A ~ O (1) admits a decomposition of its TMD quark correlator d*" 4 in terms of 
the TMD quark correlator (ft N of its nucleons, Wilson lines, and the nuclear Wigner distri¬ 
bution. The heavy nucleus in this parametric limit justifies two powerful simplifications of 
the definition (Jl]) as follows. First, the large number of nucleons justifies a mean-field de¬ 
scription of the nucleus | A(P, S )) in terms of the light-front wave functions of single-nucleon 
states \N(p, cr)). Second, because of the large number of nucleons, the initial- / final-state 
rescattering described by the gauge link U[0 , r\ is more likely to occur on the many (A — 1) 
spectator nucleons, rather than on the same nucleon from which the quark distribution is 
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taken. Schematically, we can write this as 


(A\iP p (0)U[0,r}Mr)\A) (y N -V* N ) x (N\ ^(0) u[0, r] \N) 

x (A- 1| W[0, r\\A — l) + O (kT 1/3 ) . (10) 


The resummation oi^A 1 ^ ~ O (1) 33] is a systematic way of calculating the multiple rescat¬ 
tering in the weak-coupling saturation framework, and the accuracy to leading order in A~ x ^ 
means that the sensitivity of the gauge link u[0, r] to the active nucleon is limited to O (ay). 
In this way, the TMD correlator ([I]) is factorized into a convolution of 3 factors: a wave 
function piece describing the distribution of nucleons in the nucleus, the TMD correlator of 
the nucleon itself, and a piece describing the multiple rescattering of the gauge link on the 
spectator nucleons. This essential picture should apply not only to a heavy nucleus, but 
to any system - such as a proton at high energies - in which the density of color charges 
becomes sufficiently large. 



FIG. 1. Quasi-classical factorization of semi-inclusive deep inelastic scattering on a heavy nucleus. 


For semi-inclusive deep inelastic scattering in Bjorken kinematics (see Fig.[T]), the quasi- 


classical factorization formula is 


27] 




(x, k; p, s) = Aj^—r V [ d 2+ p db~ d 2 x d 2 y d 2 k! e -dk-m*-y) W ^ P: b: S ) (11) 
(27r) 5 ^ J 

(7 

x -xp\p,o) D xy [oo~,b~}, 


where x = x^- is the quark momentum fraction with respect to the active nucleon, d 2+ p = 
d 2 p±dp + , and b M = (0 + , b~, |(x + y)) is the position of the struck nucleon. The Wigner 

































distribution of nucleons \N(p, a)) inside the nucleus \A(P, S)) is 


W a (p,b;P,S) = 


1 [ d 2 + (p- p r ) -iWb 

2(27r ) 3 J yfp + p' + 

x Y, WP, s)l Ntf, CT); X) (N(p, a); X| 4l(P, S)), 


( 12 ) 


where p = 24 p- is the average momentum of the struck nucleon in the amplitude and 
complex-conjugate amplitude and a is its spin. The Wigner distribution is normalized such 
that 


f d 2 p ± dp+ d 2 b _l db 

9+- }L / -(2^)3- P ’ 5 ) = L 


(13) 


The semi-infinite dipole scattering amplitude D xy [oo , b ] describes the final-state rescat¬ 
tering on the fraction of nucleons at depths greater than b~ and is given in covariant gauge 
by 


D Iy \oo ,b ] = — (Tr[l4[oo ,b ] Vj[oo ,b ]]), 


(14) 


where 


V x [oo ,b } = V exp 


i9 J 

. b~ 


(15) 


is a Wilson line in the fundamental representation and T a is the fundamental generator of 
SU(1V C )- The angle brackets in Eq. (TT4li denote averaging in the nuclear wave function. 

The origin of the nuclear TMD decomposition (HID is_ in a similar decomposition for 


the quark production cross section in the SID1S case 


27|, illustrated in Fig. [2] (with a 


similar cross-section analogy also valid in the DY case). There the virtual photon traverses 
the nucleus until interacting with one of the nucleons (the same nucleon in the amplitude 
and the complex conjugate amplitude). The quark produced in the interaction propagates 
through the rest of the nucleus, interacting with the nucleons. If one works in the covariant 
gauge or the A~ = 0 light-cone gauge (with the nucleus moving along the x + light cone) 
the interactions of the quark with the nucleons in the nucleus are instantaneous Coulomb 


gluon exchanges 


28] which contribute to the Wilson line describing the quark propagator. 
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FIG. 2. SIDIS cross section as a square of the scattering amplitude explicitly illustrating the 
^“-ordering of the nucleons in the nucleus. The solid vertical line denotes the final-state cut. 


The quark production cross section is 27 ] 


d a l*+ A ^q+ x 

cPkdy 


=Ag + _Y J i 

a 

r d 2 k' 

x J Wf 


dp + d 2 p db 

(2tt) 3 


d 2 xd 2 y W a (p,b]P,S) 


gfj^+N^q+X 

,-i(k-k')-(x-y) _ ( \ D [ qq — U~ 

il-ld dy P ’“ ,,[ ’ 


(16) 


Note that the light-cone Wilson lines describing the produced quark in the amplitude and 
in the complex conjugate amplitude (in the above-mentioned gauges) have to be at different 
transverse positions x±_ and y± since the transverse momentum of the quark k± is fixed (see 
e.g. (fill)- Together these Wilson lines contribute to the “staple” U in the TMD definition 
(Op). (The transverse link at infinity is zero in the covariant gauge and in the A~ = 0 
light-cone gauge.) 

Eq. (1TT1) can be simplified further by replacing the dipole scattering amplitude D xy with 
its symmetric part S xy = \{D xy + D yx ) since the anti-symmetric (odderon) part is sup¬ 


pressed by A l A 19]]. If an analytic expression is desired, this symmetric part of the dipole 
scattering amplitude can be evaluated in the quasi-classical GGM/MV multiple scattering 
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approximation to give 

S xy [co~,b~] = i (DJoo - , 6“] + £V[oo~, 6“]) 


( 17 ) 


= exp 


1 


x-y\ 2 T Q 2 ,(\‘-?\ T ) 


R-(m T )-b- 


In ■ 


4' VSVI 2 IT ' V 2 J R-(|2±K| r ) J \x-y\ T A_' 

Here Q s {bx) is the saturation scale at impact parameter bx, R~(bx) is the longitudinal radius 
of the nucleus at impact parameter b T and ( ) * s ^h e f rac ti° n of nucleons which par¬ 

ticipate in the final-state rescattering. The logarithm with infrared (IR) regulator A, which 
is often neglected, is only important for recovering the perturbative larg e-kx (small- \x^y\x) 
asymptotics. Alternatively, the Wilson line operators can be evaluated on a configuration- 


by-configuration basis if desired using Monte Carlo methods along the lines of 62 


replacing D xy with S xlJ , we can rewrite the quasi-classical factorization formula ([IT]) as 

9+ 




(x, h P, S ) = A-f-- V / d 2+ p d 2 ~b d 2 r d 2 k' _ ^)- 

(2vr) & ^ J 

x w a (p, 6 ; P, S) a) SfcX] 1 , 


63J. By 
(18) 


where we have changed variables to b = and r = x — y and shifted the integration 
variable k!_ —> kd_ + xp. 

II. 1.3. Lorentz- Covariant Spin Decompositions 


The quasi-classical factorization formula (1181) contains a sum over the spins a of the 
intermediate nucleons. In 271], a particular spin basis was chosen (the transverse |±x) 
basis) and was used to evaluate this spin sum; this effectively limited the polarizations 
of the intermediate nucleons to be either unpolarized, or polarized in the direction of the 
chosen basis. As we will now show, in a more complete treatment, the polarization of the 
intermediate nucleons can be arbitrary - unpolarized, longitudinally polarized, or polarized 
in either of the transverse directions - independent of the choice of spin basis. 

Let us first write a decomposition of the Wigner distribution (Il2li for nucleons with an 
arbitrary spin state | S') in terms of the light-cone hclicity basis |±). For simplicity, we will 
restrict ourselves to an unpolarized nucleus, but which may have polarized nucleons, which 
we write in the compact form 

W( %b,S) = J ^M,+*(y|(p) I N(p',S)} (N(p,S) I A(P)), (19) 
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where 5p = p — p ', p = ^(p + p 7 ), and (+ = (0 + ,& , b). As described in 64], an arbitrary spin 


state can be decomposed in terms of the light-cone helicity basis as 

IS) = cos 11+) +sm|e‘^|-); 


( 20 ) 


in the rest frame (R.F.) of the nucleon, the eigen-axis of spin projections for this state is 
given by the three-dimensional spin vector 


S R =' (sin 9 cos 0, sin 9 sin 0, cos 6) = (S'!, S'!, A) 


( 21 ) 


with A 2 + Sj> = 1. The Lorentz-covariant spin vector S M is obtained by defining ( 1211 ) as the 
spatial part in the rest frame and boosting it to a frame in which the nucleon momentum is 

3 


p 


Qnf \ _ ( Sr.F. ■ P X . P 

s (P) = -, S R . F . + — 

m m 


P • Sr,] 


E + m 

R = (o ,S R . F ) = ( 0 , Si, Si, A). (22) 

Note that, by construction, p fl S^(p) = 0. We have also made use of the fact that the nucleons 
do not interact with each other in the quasi-classical approximation at hand 661-168], 

By forming the outer product of ( 1201 ) and using ( 1221 ) , we obtain 


ms\ R = i 


1 cl 
2°-L 


|S)(S| R d -1 


l+X+l + l-x- 
l-X+l + l+x- 

1+) I-) 


+ ^ 


+ isi 


+ s 


l+X+l - l-x 
X+l - <I+X 
(+1 


R.F. ■ 


a 


(23) 


Eq. (l22l) constitutes the canonical definition of spin |65|,|66], obtained by boosting (0, out of the 

rest frame with a single rotationless boost. This can also be regarded as defining S' 1 = W^/m with W M 
the Pauli-Lubanski four-vector, and obtaining Sr.f. via Eq. (1221) . Alternatively, one may use light-front 
boosts, obtaining an expression for S M in terms of S' R F , which is the spin in a rest frame rotated with 


respect to the rest frame used in the second line of 


(see 


65J and the footnote 35 on page 213 of 




Ultimately, for the situations we consider (nonrelativistic nucleon motion and an ultrarelativistic boost 
along the -fz-axis), both definitions of the spin three-vector coincide. We thank Cedric Lorce for helping 
us clarify these subtleties. 
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where we have introduced the Pauli matrices and unit matrix. If we define a Lorentz- 
covariant set of Pauli matrices analogous to (122|) 


N _ / ° -V -> , P 

(?) = - ,<7+ — 

\ m m 

R.F. 


p ■ a 
E + m 


= (0,3) = (0,al,al,al) , 


(24) 


where again p p cr^(p) = 0 by construction, then we have S fl (p)a^(p) = —Sri?. ■ a in any 
frame, and 


ms\ = - L 


i+) i-) 


s i*(p) ^(p) 


(+i 


(25) 


Inserting this into (TT9|) we obtain 


W(p, b, S ) = W unp (p, b ) - s,(p) Wp 0l (p, b) 


with 


bp ?mpiPi b) — 'y ^ 


AA' 


W\y(p, b) 1 A / A 


= -T\[W(p,b)} 


(26) 




AA' L 


Wxx>(p, b) fry X (p) 


= -Tr[W(p,b) a* 


(27) 


where we have also introduced the Hermitean (2 x 2) matrix 

W xy (p,b) = J |p e - w (4(P)| N(p',\)) (N(p, A')| A(f>)>. (28) 

The projections W unp and —S p Wp 0l select out the Wigner distributions of nucleons which 
have zero polarization, or polarization S^, respectively. And although they were derived 
here in terms of the light-cone hclicity basis |±), they are invariant under a change of basis; 
any unitary rotation of the spin states transforms the matrix IT AA /, but it also rotates the 
Pauli matrices such that the traces (127)1 remain invariant. We can also write a decomposition 
of the matrix II / AA / in terms of the complete basis {!,<?} to invert Eqs. (1261) . ([27| ): 


From 


W X y{p, b) = W unp (p,b ) 


we obtain the trace identity 


Wpoi,^, b) 


J AA' 


<J M (p) 


J AA' 


--TrK(p) a"(p)]=g^- 


p fJ, p u 

m 2 


(29) 


(30) 
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which allows us to recover (|27|) from (|29|) by tracing the latter with a u or 1 and using the 
fact that P/jWph = |Tr [W p, L (7 ,l {p)} = 0. 

Note that the normalization of the matrix Wigner distribution W\y(p,b) is 


9+- 


d 2 p± dp + d 2 b± db 


■Tr[W(p,b)] — 1, 


(31) 


such that 


9+- 


d 2 p± dp + d 2 b± db 
(2vr) 3 


hhu n p(p, b) 


1 

2 


(32) 


In a similar way, the quark correlator (JTJ) of the nucleon can be expanded in a basis of spin 
states, and the polarized and unpolarized parts can be projected using the Pauli matrices: 


where 


S ) = 4> U np(x,k) - Sp^ ol (x,k) 


4* XX' (*^5 ill) fiunpfa') 

l 

4^ pol , 5 in) 

d^ 



AA' 

_ 


Punp (D —) 




0AA '(X,k) 1 A / A 
0AA '(X,k) 


]p-\[(j)(x,k)\ 
^Tr [<f>(x,k) 


(33) 


(34) 


and the matrix form is 

*«,(*,£) = J±=j Jte-re^ (JV(P,A)| ^(0)M[0,r]V> o (r) |JV(P,A')). (35) 

Note that Sp and d M above depend on the (averaged) momentum of the nucleon p, which 
is not shown explicitly. The expansion (J3]) of the quark correlator into the 8 leading-twist 
TMDs is defined in terms of the longitudinal spin Sl and the transverse spin S_. The Dirac 
projections of the correlator 0^ as in (|9]) are boost-invariant, so we can evaluate them in 
the rest frame in which S = (S(J_, Sj_, Sl)', this allows us to identify the projections (pulp and 
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^pol ] i 11 terms of the TMDs as 


= /i 


unp 


C R = o 
Cf *= (» 


-^+_Tr[ 7 r] 


Li 

+ [^hi 

m 


-c; + _Tr[7 ±i 7 T] 


:9-\ —Tr[ 7 5 7 _ r] 


+ 


k\ 


±h i L 


m 


r3+-Tr[7 5 7±*7 r] 


, i [T] JLF. 

y poZ, ± — 


Li 

ft -L ,U f-L 
fc T J IT 
m 


-^+-Tr[ 7 r] 


,ki 

+ 9 IT 

m 


- 9 + _Tr[ 7 5 7 -r] 


Li L3 

8 lj h lT + -±-±h 


m 


2 "IT 


r 9 +-Tr[ 7 5 7±i7 r] 


(36) 


In the derivation of (1 1 81) . a sum over the spins of the intermediate nucleons occurs inde¬ 
pendently in the amplitude and complex-conjugate amplitude (see Fig. Q]). Thus the spins 
entering the quasi-classical factorization formula are necessarily non-diagonal, which can be 
compactly expressed in terms of the matrices Ww and (j)\ 


$ A {x,k P) = A 


= A 


9+- 

(2vr) 5 


^ f d 2+ pd 2 bd 2 rd 2 k'( 

XX' J 


i(k—k'—xp)-r 


X W^(j>,ly,P)^, x (x,k'-,p)Sp T 2] ] 
9+ ~ [ d 2+ pd 2 ~b d 2 r d 2 k' e -*(£-£-*p)-r 


(2vr 


x Tr 


W(p, b ; P) <f>(x,k';p) 


o[oo ,b ] 

(r T ,b T ) ’ 


where we have restricted the discussion to an unpolarized nucleus. From 
(USD, the trace of the two matrices gives 

1 


(37) 

and 


-Tr 


W ■ (j) 


1 Vunptftunp fFpo/. // tppol 


R.F. 


1 f j/np tftunp + bFpoZ 


(38) 


which sums over all four independent spin states of the intermediate nucleons. Thus we 
recover a form similar to (fl8jl in which the spins of the nucleons are diagonal (i.e., the same 
in the amplitude and complex-conjugate amplitude), 

$ A (x, k- P) = 2 ' 4 ^7 [ d 2+ p d 2 ~b d 2 r d 2 k' 

(27t) & J 

x (w unp (p, 6; P ) <j> unp {x,ti\p) - W po i,p(p, b; P)(f>p 0 i(x,k!]p)) s1 ™ t / t) \ (39) 
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as illustrated in Fig. [3j The sum now effectively runs over the four polarizations: U = 
unpolarized, L = longitudinally-polarized, and T- 7 = transversely polarized in the j direction. 
This formulation makes the projections defined by (1251) . (1551) independent of the choice of 
spin basis and expresses the spin dependence in a manifestly Lorentz-covariant form. 



FIG. 3. Schematic decomposition of the spin sum occurring in (1391) . Summing over the 4 indepen¬ 
dent components of the (2 x 2) matrices [</>]aa'> [W]a'a corresponds to summing over the 4 possible 
intermediate polarizations of the nucleons: U = unpolarized, L = longitudinally polarized, T J = 
transversely polarized in the j direction. 


II.2. Parameterization of the Wigner Distribution 

The quasi-classical factorization formula (1591) describes how the properties of the Wigner 
distribution and the multiple rescattering build up the TMDs of the nucleus from the various 
TMDs of the nucleons. One powerful feature of (|39l) is that the Wigner distribution is con¬ 
structed purely from the light-front wave functions (see, e.g. (125]) ). without contamination 
from the scattering dynamics embodied in the gauge link. It is therefore invariant under 
several symmetries which the TMDs themselves are not, such as PT. We would like to use 
all the applicable symmetries to parameterize the types of structures which can occur in the 
Wigner distribution. 

The Wigner distribution also does not know about the direction of the collision axis; 
that information enters operatorially through the gauge link 5/[0,r] or diagrammatically 
through the interaction with the virtual photon. We may choose to describe the functional 
dependence of the Wigner distribution using variables which are appropriate for a collision 
along the £-axis (e.g. p + ,p), but the distribution itself should not treat z as a “special” 
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direction. Naively, one would expect that, when viewed from the rest frame of the nucleus, 
the Wigner distribution W(p, b ) should possess three -dimensional rotation symmetry due to 
the lack of a preferred collision axis. This is in contrast to the factors 4>(x,k) and S(tt , hr) 
containing the gauge link, which have at most two-dimensional rotational invariance due to 
the special role of the longitudinal direction. 

There are some subtleties about the application of full rotational symmetry to the Wigner 
distribution constructed from light-front wave functions; these are addressed most clearly 
using the “covariant light-front formalism” of 69]. We discuss these issues in Appendix [A] 
The conclusion is that, if the nucleons move nonrelativistically in the nucleus, then their 
wave functions - and hence the Wigner distribution in (1A19I) and (IA20D - do possess the 
desired rotational invariance in the nuclear rest frame. The boost-friendly variables ct, b~ 
are given in terms of the rotation-friendly variables p, b through (IA21D . We can now take 
advantage of this symmetry to constrain the functional form of the Wigner distribution. 


II. 2.1. Wigner Distribution of an Unpolarized Nucleus 


For an unpolarized nucleus with a non-relativistic distribution of nucleons, the Wigner 
distribution (jA20j) possesses manifest rotational invariance in the rest frame and depends 
on the nucleon spin only linearly. Moreover, because the Wigner distribution is built purely 
from the nucleon wave functions, it possesses the unbroken discrete symmetries of (ordinary) 
parity and time reversal. The P and T symmetries of the wave functions translate into the 
Wigner distribution in the expected way: 

W(p, b, S) £ w(-p, -b, 8) = w(-p, b, -S). (40) 


We can therefore constrain the form of the Wigner distribution using these discrete symme¬ 
tries and rotational invariance. 

For the unpolarized distribution, requiring rotational invariance yields 

p 2 ,b 2 .p- b 


W unp (p,b) = W ., 


unp 


(41) 


but the quantity p ■ b is odd under time reversal, so W unp can only depend on it through 
even powers: 


W unp (p, b) = W unp ]?,b 2 ,(p-b) 


(42) 
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Now let us change variables to the boost-invariant quantities using (1 A21 I) : 


f=p 2 T +(a- \) 2 M 2 a 

b 2 = b 2 T + (^) 2 (g^P + b-r 

(. P-b) 2 = [p-b-(a- \)(g+„P + b-)] 2 

= (P • kf T + (a - \) 2 (g + -P + b-) 2 - 2(p ■ 6)(a - i)(<? + -P + r), (43) 

where a = p + /P + is the longitudinal momentum fraction of the nucleon in the nucleus. In 
terms of these quantities, the unpolarized distribution depends on 

W unp (a,p; b~, b) = W unp [p 2 T , b 2 T , (p ■ bf T ; (a - \) 2 , ( P + b~) 2 ; (p • 6)(a - i)(P + r)] , (44) 

where we have separated out the dependence on the longitudinal and transverse variables. 

There is one interesting term which appears to mix the longitudinal and transverse vari¬ 
ables, arising from the cross term (p\\b\\)(p-b) of ( p-b ) 2 . In principle such a term is not prohib¬ 
ited by the symmetries of Wigner distribution; however, looking back at the quasi-classical 
factorization formula fl39]h we see that this term will not contribute to the cross-section. The 
reason is that the only other dependence on b arises from the dipole scattering amplitude 
' (TTTT) . This factor does not possess the same three-dimensional rotation invariance as 
the Wigner distribution (1A20I) . because the multiple scattering depends on the depth b~ of 
the struck nucleon. It does, however, retain a residual two-dimensional rotation invariance 
in the transverse plane, depending only on b~, b Thus a term like J d 2 b (p\\b\\)(p ■ b) fibx) 
would integrate out to zero in (|39|) because of antisymmetry under b —> —b. Therefore the 
terms of the Wigner distribution which survive the d 2 b integral can only depend on the 
square of this cross-term, (p ■ b)^(a — l/A) 2 (P + b~) 2 } and the dependence on these squared 
factors are already taken into account. Thus we can replace 

W unp (a,p-,b~,b) W unp [pT,b?r, ( p-b )|; (a - \) 2 , ( P + b ~) 2 ] (45) 

without missing out on any terms which would survive the d 2 b integral in ()39lh We use an 
arrow here rather than an equality to emphasize that we are now dropping contributions 
to the Wigner distribution which are permitted by symmetry, but would not survive to 
contribute to the nuclear TMD. 
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In fact, we can carry this argument one step farther: after the integration 


I d 2 bW unp (a,p]b ,b)S(r T ,b T ;b ) 

is performed, all sensitivity to the direction of the two-vector b drops out. Thus we can 
replace 

f-i<4 -> (46) 

in the Wigner distribution now, without loss of generality in the types of terms which can 
contribute to the quasi-classical factorization formula (l39]h Doing so further reduces the 
number of terms which can contribute, since 

ip ■ —)t = (&!l&±) ->■ \pr h T- (47) 

Thus, from all of these simplifications, we can write the relevant part of the unpolarized 
Wigner distribution as 

W unp (a,p;b-,b ) =► W unp [p 2 T ,b 2 T ; (a - ^) 2 , ( P + b~ ) 2 ] . (48) 


By considering the constraints due to parity, time reversal, and three-dimensional rotation 
invariance of the Wigner distribution, together with the two-dimensional rotational invari¬ 
ance of the multiple scattering factor (gauge link) to which it couples, we have reduced the 
unpolarized distribution down to a form which is manifestly even under p —y —p, b —> —b, 
(a — jf) —> —(a — 4-), and b~ —> —b~. These symmetries strongly constrain the interplay 
between the factors W(p, b), 0(x, k'), S(r, b) entering the quasi-classical factorization formula 

m- 

Similarly, we can impose rotational invariance on the polarized part of the Wigner dis¬ 
tribution: 


S-W pol (p,b)= (5-6) W 1 [fP,P,p-b\ + ( 

s-p^j w 2 

f,b 2 ,p- b 

+ (s-(6xp)) W 3 

f,b 2 ,p- b 



(49) 


The spin-dependent factors (S-b) and (5 -p) are odd under parity, while all of the arguments 
f,b 2 , p-b are P-even; therefore these types of spin dependence cannot enter into the polarized 
Wigner distribution. The only spin dependence which is permitted by parity is the factor 
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(S ■ L) describing spin-orbit coupling, where L = b x p is the orbital angular momentum of 
the nucleon. Since (S ■ L ) is even under both parity and time reversal, the T-odd quantity 
(p ■ b ) can only occur in even powers, just like in the unpolarized distribution: 


S-W pol (p,b) 


S-(bxp))W 3 f,b 2 ,(p-bf 


(50) 


Now let us again change back to the boost-invariant variables ( 1A21jl . 
expressions in (T43j) . we must account for the spin-dependent factor 


In addition to the 


S • (6 x p) = S\\{b x p) + b\\(p xS) + p\\(S x b) 

= A (6 x p) - 7+ -P + b-) (p x S) + (a - \)M A (S x 6) (51) 


to write 

-S A W^(a,p-,b-,b) = (A (6 x p) - jj^(g+-P + b~) (pxS) + (a- j)M A (S x 6)) 

x W 3 [p 2 T ,b 2 T ,(p-b) 2 T ; (a-}) 2 ,(P + b~) 2 ; (p- b)(a - i)(P + b~)] . 

(52) 

As before, we now replace b l ± tP ± —* \b^8 1 ^ to keep only the terms which can contribute to 
the nuclear TMD. In addition to reducing (p ■ b)^ —> ^Pt^t-, this also simplifies the spin 
dependence for the factors which couple S_ to b: 


[A (6 x p)] [(p ■ 6)(a - \){P + b )] ->■ |A(a - \){P + b ) (p x p) = 0 
[(« “ i) (S x k)] [(p ■ b)(a — \)(P + b~)] ->• -\(a - \) 2 ( P + b~)(p x S). (53) 

We see that the spin dependence of the form A (bxp) has dropped out completely, and the 
spin dependence of the form {a — \)(S_x b) has reduced down to the form ( P + b~)(p x S_). 
Thus, without loss of generality, we can identify the structure ( P + b~)(p x S) as the only 
spin dependence which survives the d 2 b integral to contribute to the nuclear TMDs: 

-S„W^(a,p'b-,b) =c (p x S) W 0 AM[p 2 T,b 2 T ; (a - A) 2 , («-)’] . (54) 

Again, consideration of the relevant symmetries has reduced the form of the polarized Wigner 
distribution significantly. The structure that survives contains a prefactor linear in the 
transverse spin which is odd under b~ —y —b~ and p —> —p, times a function which is 
manifestly even under p —$■ —p, b —> —b, (a — —> —(« — -j), and b~ —> —b~. The 
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FIG. 4. Spin structure of the unpolarized nucleus. An unpolarized nucleus can give rise to unpolar¬ 
ized nucleons or transversely-polarized nucleons, which contribute to the quark distribution through 
various nucleonic TMDs. Interestingly there is no contribution from longitudinally-polarized nu¬ 
cleons in the unpolarized nucleus (1561) . 


designation “OAM” reflects the fact that this structure originated from the presence of L ■ S 
coupling in the rest frame. 

Altogether, (I48p and (I54p allow us to parameterize the most general structure of the 
Wigner distribution which is consistent with parity, time reversal, and 2D and 3D rotation 
symmetries. The structure which can contribute to the quark TMD of an unpolarized 
nucleus is 


W(a,p;b~,b; A,5) =4- W„ p [^.,4; (a - i) 2 , (P + !T) 2 ] 

+ (tt) (£ x 3 W oam [pt> b 2 T- (a - \) 2 , (P + fT) 2 ] , (55) 

or, in the notation of (l26lh 

w uv (p, b) = W„ p [p 2 T ,b 2 T ; (a - i) 2 , (P+fT ) 2 ] 

W£,(p,6) 0 

R = 0 

.(p- 4 ) h = ( i± =Sr~)'P'i. t T W OAM[PT,^T^ (a- y 2 .(P + 6”) 2 ] ■ (56) 

Eq. (1551) is illustrated in Fig. [4} The parameterization given here was derived for the distribu¬ 
tion of nucleons with polarization Sjy = (S_, A) moving nonrelativistically in an unpolarized 
nucleus, for which the only nontrivial spin dependence entered as (L • Sn). It is straight¬ 
forward to extend this parameterization to describe a polarized nucleus as well. In addition 
to the equivalent spin-orbit coupling (L • Sa ) which was the source of the “orbital angular 
momentum channel” to generate the nuclear Sivers function in j^ . there are spin-spin 
couplings like (Sn ■ Sa) and spin-orbit-spin couplings like (Sn ■ p)(p • Sa)- We will leave 
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this generalization to the full structure of the Wigner distribution of a polarized nucleus for 
future work. 


II.3. Quasi-Classical TMDs of an Unpolarized Nucleus 

The TMD quark correlator & A np of an unpolarized nucleus from (foil contains just two 
leading-twist TMDs: the unpolarized quark distribution f A and the Boer-Mulders distribu¬ 
tion hi A : 


= (fi{x,k T ) 


-3 + _Tr[7 T] 


+ 


k\ 


± h^ A (x, k T ) 


M a 


r ^+_Tr[7 ±i 7 T] 


(57) 


Because f A corresponds to unpolarized quarks and h^ A corresponds to transversely-polarized 
quarks, these two TMDs can be projected from the correlator as in fl5]) by tracing the Dirac 
structure with T = 7 + or T = 7 5 7 + 7^_, respectively: 


hr) 


k\ 


ji "'-L u±A 


M a 


hf A (x,k T ) 




(58) 


These nuclear TMDs are related to the TMDs of the nucleons through the quasi-classical 
factorization formula 




2 A g + - 
(27t) 5 


d 2+ p d 2 ~b d 2 r d 2 k' eMk-H-ipyr 


X ( w unp (p , b ) Cip(£, U) - W po i^{p, b)4>^ ] (x, U)) S [ ™ T / T ) ] , (59) 


with the Wigner distribution of nucleons parameterized by (1561) . Since the nucleons can be 
polarized, the active quarks can come from one of several nucleonic TMDs depending on the 
polarization, as in (136|) . Combining the Wigner distribution for a given nucleon polarization 
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with the associated TMDs, we obtain 


W unp (p, b ) <i>N p {x, U) = W unp [p 2 T , b 2 T ; (a - \) 2 , ( P + b ) 2 ] 

'1 


x i ( A (A ^t) 


: 9 +- Tr[7 A 


+ 


h' 1 

K ± u±N 


m N 


hf (x, k' T ) 


:9+- Tr[7±*7 r ] 


-) pie^ Woam [p't, bl; (a - A) 2 , (.P+fr ) 2 ] 

u'i \ ri 

ft/ ii , I M / / , / x \ -L 


_ f g+^P+fe- 
AI a 


X 


■”-L ,*7 f-LAf | 
t TJlT 1 

m N 


x, k' T ) 


4 


^+_Tr[7 T] 


+ 


k' j 

—g^(x,k' T ) 

m N 


-(? + _Tr[7 5 7-r] 


AAf T (f, *£,) + 


U'iU'J 

^Lh^(x,k' T 

m N 


-^+_Tr[7 5 7 ±i 7 T] 


(60) 


where a = p + /P + and x = x/a. The projection [T] of the nuclear TMD quark distribution 
<& A in (l58l) selects the distribution of quarks with a given polarization ( U,L,T J ). This 
projection also applies to the nucleonic TMD quark distributions </> in (1601) . picking out a 
specific combination of nucleonic TMDs which can contribute to a given TMD of the nucleus. 
The orbital dependence on the quark momentum A/ carried by each of the TMDs can couple 
to the orbital dependence on the nucleon momentum p in the Wigner distribution, giving rise 
to nontrivial interplay between the TMDs of the nucleons, the orbital angular momentum 
these nucleons carry within the nucleus, and the resulting TMDs of the nucleus. Now let us 
use these properties to separately study the quasi-classical decomposition of the unpolarized 
quark distribution / A and the Boer-Mulders distribution h A A of the nucleus. 


II. 3.1. The Unpolarized Quark Distribution f A 

The projection T = 7 + in (158]) selects out the unpolarized quark distribution ff in the 
nucleus and the unpolarized quark distribution /A and quark Sivers function /A V of the 
nucleon in (1601) (see also (j9]) ) . Using this in (l59j) gives 



x ( 'w mr (p, b ) P(x, k’ T ) - Tg^(P + b-)(p ■ k') Woam(p, b) k’r) j , ( 61 ) 

which shows that (L- S N ) spin-orbit coupling given by Woam can result in a mixing between 
the unpolarized quark distribution ff of the nucleus and the quark Sivers function of 
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the nucleons. This is the mirror effect to the OAM channel that was found in [27[, in which 
(L ■ Sa ) coupling in a polarized nucleus gave rise to a mixing between the Sivers function 
/fr 4 of the nucleus and the unpolarized ff distribution of the nucleons. 

The unpolarized quark distribution f\ is a PT-even function and therefore process in¬ 
dependent (see [70;] and others); that is, it should give the same result when calculated 
with a future-pointing gauge link like (ITTj) appropriate for SIDIS as when calculated with a 
past-pointing gauge link like 


q[b , —OO ] 

(rr,br) 


exp 



( b- +fl-(6r) \ 

V 2 R-(br) ) 



(62) 


appropriate for the Drell-Yan process (DY). In the quasi-classical approximation, the dif¬ 
ference between the future-pointing gauge link of (ITT)) and the past-pointing gauge link of 
(l62|) is in the fraction of nucleons which contribute to the scattering. In the future-pointing 
case, the nucleons at depths x~ greater than b ~, ( b~ < x~ < R~{b T )) contribute to the 
final-state rescattering; in the past-pointing case, the nucleons at depths x~ less than b~, 
(— R~(br) < x~ < b~) contribute. For a uniform distribution of nucleons, this is what gives 
rise to the factors 


R-(br) 


T(b 7 


dx pn(x , br) = 


R~(bi 


T(b 7 


dx Pn{x , b T ) = 


-R-(br) 


2 R~(br) 

b- + R-(br) 
2R~(b T ) 


(63) 


describing the fraction of nucleons which participate in the scattering. Here T(br) is the 
number density of nucleons per unit transverse area at impact parameter bx- Because ff 1 is 
PT-even and process-independent, it should give the same result whether evaluated with the 
future-pointing gauge link of (fTTj) or the past-pointing gauge link of ()62|) . This is clear for 
the trivial channel because both ff and ff are invariant under the PT transformation, 
and the longitudinal b~ integral is also preserved 


db e 


_1.2/i2/ R ( b T )-b 
4 rVsl 2 R~(b T ) 


In 


rrA w unp [(P + b-) 2 ] = (b- (-6-) 


= / dP 


-\rfQl(tPPLJtfR | in-i 


2R ~ (bT) ) rrA W unp [(P + b~) 2 ] (64) 
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by a simple change of variables. The Sivers function /f^ v , on the other hand, is PT-odd 
and changes sign between the two processes, but its associated longitudinal b~ integral 


db e 


_1_2 n 2 ( R (b T )-b 
4 r TVsl 2 R-(b T ) 


In 


VtA ( P + b ~) Woam [( P + b -) 2 ] = (b~ (~b~) 


_1 r 2 Q 2('h+TM'l 

db e 


In ■ 


-(M ; r T A (p+r) [(p+ 6 -)2] (65) 

also changes sign. Thus both channels of (j6Tj) are overall even under PT, yielding the 
process-independent nuclear TMD as required. Similar to [27], the (L ■ Sn) spin- 
orbit channel here requires at least one additional rescattering to be nonzero; in the limit 
r^Q 2 1, corresponding to either large transverse momentum or low densities, the multiple 
scattering factor S(rr, hr) can be neglected, and the b~ integral (Rj 51) vanishes. The depth 
dependence on b~ of the spin-orbit factor ( P + b~)(p x5), coupled to the angular dependence 
of the Sivers function, provides a second PT-odd factor (P + b~)(jp ■ }z) which enables a PT- 
odd nucleonic TMD like to contribute to a PT-even nuclear TMD like f A . As we will 
show, the same (L ■ Sn) spin-orbit coupling in the nucleus also generates contributions to 
the PT-odd nuclear Boer-Mulders function from various T-even nucleonic TMDs. 

To better understand our result, let us evaluate Eq. (jfiTl) for specific models of the Wigner 
distribution. First let us consider the case when the internal motion of the nucleons in the 
nucleus is negligible, that is the nucleons are static with respect to the nuclear center-of-mass 
(the “standard” GGM/MV model in saturation physics). The corresponding unpolarized 
Wigner distribution is (cf. Eq. (65) in [ 2 ^) 


W unp (p,b) R = ~ 6 2 )6' 3 (p) 


3n 


fGmp(p, b) 


R 3 

9 + ~ (2 7r) ; 

2 A 


p(b, b )S (p)8[p 


P + 

~A 


( 66 ) 


(67) 


with the constant nucleon number density (per d 2 b± db ) 

61 ) na. 

(The normalization for W unp is defined in Eq. (1H2]) .) Since the nucleons are static, Woam — 0 
and the second term on the right-hand side of Eq. (16T1) vanishes. Finally we take for the 
unpolarized quark TMD of the nucleon the lowest-order perturbative expression for a quark 
target 48, 57] 

a s CV 1 + x 2 


fi hr) = 


2 7T 2 kj, 1 — X 


( 68 ) 
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which we for simplicity consider in the kx limit. As usual Cf = {N 2 — 1)/(2N C ) is 

the fundamental Casimir operator of SU(A r c ). Using ([68]) . (166|) and (ITT)) in Eq. (16T|) yields 


/i 4 (u hr) = 


N c 1 + x 2 


d 2 bd 2 rea¬ 


cts 47r 4 1 — x J 
The low-rc (and low x = Ax) limit of Eq. 


r T L 


exp 


■J 1uT_ 


(69) 


calculation of a similar quantity in 
unintegrated quark distribution in 


is in complete agreement with the previous 


H (see Eqs. (22) and (23) there keeping in mind that 
is defined per dkpdx/x unit of phase space, while 
TMDs are conventionally defined per d 2 kxdx). The quark saturation scale in Eq. (I69|) 
is Q 2 (bx) = / l / na 2 T(bx) Cf/N c . Note that the agreement between our formula (16T|) for 
the unpolarized quark distribution with the results obtained from a different method in 
48] validates our approach and suggests that other observables, like the unpolarized gluon 
distribution, when calculated in our technique would agree with the results in the literature 
|l8, 32, 61, 711. 

As is well-known, putting the logarithm in the exponent of Eq. (I69|) to one (as a negligibly 
slow varying function) casts the equation in a form where the integral over r >r can be carried 
out analytically. This gives 

N c 1 + x 2 


A (x, hr) = 


d b T 0, 


fa? 

qWi 


(70) 


a s 47t 3 1 — x 

This approximation is valid only for kx not much larger than Q s . 

To get a feeling for the correction to Eq. (1(7)|) resulting from the spin-orbit coupling term 
in Eq. (1611) assume that 


W 0 AM(p,b ) fa £ 


9 + (2 7r) ; 

2 A 


P(k,b ) 


u-PtA 


7r m 


N 


5 \P + -- 


P+ 

A 


(71) 


where £ is a parameter responsible for the strengths of spin-orbit couplin g an d we neglect the 
small change in p + due to the nonrelativistic orbital motion. Just like in 27] we approximate 
the nucleon Sivers function by the quark target expression, assuming again that kx rriN'- 
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Employing Eqs. (ITT)) . ([72]) and (fTTl) in Eq. (I6TT) we get the following correction to the unpo¬ 
larized quark TMD of a large nucleus: 
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where p = is the nucleon number density in the rest frame and we have made use of 
the -C Q s condition and put ln(l/r T A) ~ 1 in the exponents. Comparing Eqs. (T73|) and 
dZQD we see that the spin-orbit contribution is parametrically different by a factor of £,/a s . 
Without the knowledge of the (most likely) non-perturbative parameter £ it is impossible to 
say whether the factor of £/a s indicates suppression or enhancement. If we take £ ~ a 4 by 
analogy to atomic spin-orbit coupling in quantum electrodynamics (QED), then £/ct s ~ 
which is parametrically very small. However there is no good reason to justify such a choice 
in QCD. Further work is needed to fully quantify the role of spin-orbit coupling in the 
unpolarized quark TMD. 


II.3.2. The Boer-Mulders Distribution 


±A 


Similarly, the projection T = 7 5 7+7^_ in ([58]) selects the Boer-Mulders distribution hj~ A 
in the nucleus and the functions hj ~ N , h ^ T , 44 of the nucleon in (l60|) (see also ([9])). Using 
this in (I59|) gives 
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It is convenient to separate the diagonal part 8^ of the quantity in parentheses, 

u'o yt y2 / yj ye i \ 

8 je h4(x, k' T ) + ~ L 2 ~ L hjj v (x, k ' T ) = 8 jl 4 (x, k' T ) -\ -1- ( + + - -8 je ) 4t( x, k ' T ), (75) 
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is the transversity TMD and 44 is the pretzelosity. Then 
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Then we can contract the free index j on both sides with — and solve for the nuclear 
Boer-Mulders function to obtain 
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(78) 


Thus we see that the same L ■ Sn coupling which allows the quark Sivers function 
of the nucleons to mix into the unpolarized quark distribution /j 4 of the nucleus in (16T1) 
also allows the transversity hf and pretzelosity of the nucleons to mix into the nuclear 
Boer-Mulders function h^ A of the nucleus in (l78]h The Boer-Mulders function, like the 
Sivers function, is a PT-odd function that changes sign between the future-pointing gauge 
link (1X71) and the past-pointing gauge link fl62j) . The trivial channel which builds hf A from 
hf N already has the necessary properties under PT, and the longitudinal b~ integral is the 
same as in (J64]). But there are also PT-even channels - the nucleonic transversity h± and 
pretzelosity - which mix into the PT-odd nuclear Boer-Mulders function through the 
role of orbital angular momentum. Just like in the case of the nucleonic Sivers function 
/mixing into the unpolarized quark distribution f A of the nucleus, the necessary reversal 
of PT symmetry is provided by the depth dependence ( P + b ~) in the longitudinal integral 
(]65]h And just like before, this mixing requires at least one rescattering to provide this extra 
PT-odd factor, so it vanishes in the limit of large kx or low charge densities. 


III. EVOLUTION 

Our goal in this Section is to include the quantum evolution corrections to the quasi- 
classical TMDs like those calculated above. We will separately consider the cases of x = 0( 1) 
(large-x) and x 1 (small-re). 


III.l. Large-x Evolution 


The TMDs calculated above and in 



were found for x ~ 0(1) \ here we need to find 


the quantum evolution corrections to a generic TMD given by Eq. (I39]h Note that we are 














working in the s ~ Q 2 k\ regime, where s is the center-of-mass energy squared per 
nucleon; hence our x in this section (and above) is order-one, but does not closely approach 
unity. Let us begin by concentrating on summation of leading logarithms of energy, that 
is on powers of cu, In s/k ^). To do this we can use the established formalism of saturation 


physics |40|-l47|,l72|-l7^. 


We work in the frame where the incoming nucleus has a large P + momentum, P M = 
(P + , ^g + ~-pr, 0), while the incoming virtual photon has a large positive q~ along with a 
comparable but negative q + momentum, = (— q~, 0). (For definitiveness we 
consider the SID1S process, though all our discussion and conclusions apply to DY as well.) 
We will work in the A~ = 0 light-cone gauge. 



FIG. 5. An example of quantum evolution corrections to the SIDIS process and for corresponding 
quark TMDs. 

The typical quantum evolution corrections to a SIDIS amplitude are shown in Fig. [5j 
Our choice of A~ = 0 gauge prevents gluons from being emitted or absorbed by the partons 
in the incoming nucleons or by the struck quark (carrying momentum k! — q in Fig. [5]); 
all these particles move predominantly in the P + -direction and do not emit gluons in the 
A~ = 0 gauge in the leading-logarithmic approximation. (We have also checked by an 
explicit calculation that the emissions off of the k' — q quark do not generate logarithms 
of energy.) Therefore, gluon emission and absorption is limited to the outgoing quark. 
Additionally, logarithms of energy can be generated (in A - = 0 gauge) only by emissions 
which happen over long periods of time; therefore, gluon emissions inside the nucleus do not 
lead to logarithms of energy. We are left only with the gluon emissions and absorptions by 
the outgoing quark after it exits the nucleus, as shown in Fig. [5j 

Such emissions are easy to sum up. Using crossing symmetry, the correction to a light- 
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cone Wilson line at x±_ in the amplitude and another light-cone Wilson line at y± in the 
complex conjugate amplitude, as appears in the cross section, can be accounted for by 
calculating corrections to a pair of Wilson lines in the amplitude: one at x± another one at 
y |_. This is illustrated in Fig. [Sand is the basis for expressions like Eq. (13Uj) . though in the 
quasi-classical case one resums multiple rescattering diagrams. The crossing symmetry in 
Fig. [6] is valid both for multiple rescatterings and for gluon emission 75]. 







FIG. 6. An illustration of the crossing symmetry for the quantum evolution corrections on a pair 
of semi-infinite Wilson lines. The thick vertical band indicates the nucleus (the shock wave) with 
all the multiple scatterings on its nucleons. 


We conclude that all we have to do is evolve a semi-infinite dipole stretching in the x~ 
direction from 0 to +oo. (In the DY case the dipole would stretch from — oo to 0, but the 
evolution equation would be the same.) In the leading-logarithm of energy approximation 


the evolution for this object is simply given 


small -x evolution of an infinite dipole 40 


by the 
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half of) virtual corrections to the standard 
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(x — z ) 2 (z — y) 2 

with the initial condition, S^oo - , b~](Y = 0), given by Eq. ((171) above. Here d Y = d/dY. 
The rapidity variable is given by the logarithm of energy, Y = In[s (x — y) 2 ], where the IR 
cutoff (x — y) 2 is beyond the control of the leading logarithmic approximation considered 
here. All we know is that this logarithm of energy has a cutoff of the order of fc T . 

Eq. d79]) is easy to solve 72]. Integrating over z± on the right-hand side yields 

(x - yf 


d Y S xy [ oo ,b ](Y) = Cf In- 


7T 




S X y[ oo ,b ](Y) 


(80) 


with p some ultraviolet (UV) cutoff having dimensions of distance. Since the shortest trans¬ 
verse distance in the problem is of the order ofl/Q~l/s (since we assume that s ~ Q 2 ), 
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we make another approximation: keeping only the leading logarithms of Q 2 we replace 
p —> 1/Q. In addition, remembering that we have s ~ Q 2 we see that 


Y = ln[s (x - y) 2 } = In — + ln[Q 2 (x - y ) 2 ] 


s~Q 2 


ln[g 2 ( x-y ) 2 ]. 


(81) 


In the resulting double-logarithmic approximation (DLA) (in InQ 2 ) we rewrite Eq. (1801) as 
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dp 2 a s Cp 


V Qo 




71 


In [p 2 ( x-yf 


S xy [oo-,b-}(Ql) (83) 




with the initial conditions S' xy [oo ,b ](Qq) Qo = 


a 


given by Eq. ffTTjl . The exponent 


39], previously derived in the saturation 


■fer vh 


in Eq. (1831) is the well-known Sudakov form-factor 
literature in 76979], For gluon TMDs, where the Wilson lines are adjoint, one would have 
to replace Cp —$■ N c in Eq. (l83j) . While Eq. (183|) is written for a fixed coupling constant a si 


running coupling corrections to it can be included following 


80 


Note that while we started out this derivation by studying evolution of a semi-infinite 
dipole with energy, we have been working at large-a: where s ~ Q 2 and evolution with 
In(s/fcy) is equivalent to evolution in In (Q 2 /kp), as follows from Eq. (j8TT) with the transverse 
momentum of quarks in the TMDs being approximately kp ~ l/|x —?/|. The Sudakov form- 
factor (183|) that we have obtained is a result of Q 2 evolution of the TMDs at large-x, and 


is different from the well-known Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
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85 ] (Revolution of the parton distribution functions (PDFs). The difference is ultimately 
due to TMDs and PDFs being different objects, with kp of the TMDs integrated over to 
obtain PDFs: in [^6| DGLAP evolution has been re-derived for PDFs defined as kp -integrals 
of TMDs. 

To summarize, the QCD-evolved large-a; quark TMDs of a large unpolarized nucleus in 
the saturation picture can be obtained from the quark-quark correlator 
2 Ag + . 
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with the semi-infinite dipole scattering amplitude given by Eq. (j83]) in the DLA and 
S xy [oo~,b^](Ql) with Ql = given by Eq. (fT71h For gluon TMDs the Sudakov form- 

factor would only be different by the Cf —> N c substitution. 


III.2. Small-x Evolution 

III. 2.1. Evolution of Unpolarized-Target TMDs 

Now let us consider TMD evolution in the small- x regime, s >Q 2 > k We begin with 
the unpolarized proton or nucleus TMDs, concentrating on the unpolarized quark TMD 
/-j 4 first. As usual, the dominant contribution at small-x comes from diagrams which are 
order-ct s suppressed compared to the channel shown in Fig. [2] above. Concentrating on the 
SIDIS process again, the dominant small-a; contribution is pictured in Fig. [71 Comparing 
Fig. [7] to the lowest-order (no multiple rescatterings) part of Fig. [2] one immediately sees 
that the former is order-a s suppressed: this is why the channel in Fig. [7] was not considered 
above. However, at snmll-x the channel in Fig. [7] is dominant, being enhanced by a power 
of 1/x compared to the channel in Fig. [2j 



FIG. 7. A leading-order contribution to the unpolarized-target SIDIS process at small-x. The 
tagged quark is marked by a cross. 

One may ask why the diagram in Fig. [7] was not included along with the evolution 
corrections considered in Sec. IIII.ll Indeed the quark loop in Fig. [7] generates a factor of 
In (Q 2 /k^), such that the suppression of Fig. [7] compared to the lowest-order part of Fig. [2] 
is diminished to a factor of a s In (Q 2 /k^). However it is well-known that the quark loop 
cannot generate a logarithm of energy; hence the diagram in Fig. [7] is a single-logarithmic 
correction, which is beyond the DLA accuracy of the evolution in Eq. (1831) and can be 
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neglected at large-x. 

The splitting of a virtual photon into a qq pair in Fig. [7] happens typically long before 
the interaction with the target nucleus, denoted by a shaded oval in Fig. [TJ Moreover, 
in the unpolarized SIDIS case all interaction with the target is eikonal and, hence, spin- 
independent. No interaction with one given nucleon is a special “knockout” interaction, 
where spin-dependence could be transferred from the target to the probe, as in [27]. There¬ 
fore the unpolarized SIDIS cross section at sniall-x, and the corresponding TMDs, do not 
have the b~ dependent effects akin to those at large-x as shown in Eq. (139|) leading to TMD 
mixing in Eqs. (1BTT) and (1751) . We also do not need to worry too much about the details of 
nucleon dynamics encoded in the nuclear Wigner distribution. 


\ 



\ 



\ 


FIG. 8. Diagrams representing the sum of all-order contributions to the unpolarized-target SIDIS 
process at small-x. The tagged quark is marked by a cross, the thick vertical band represents the 
target nucleus (the shock wave). 


With this in mind, summing up the diagrams in Fig. El where the thick vertical band 
represents the shock wave and we explicitly show the four contributions where the 7 * —>■ qq 
splitting happens either before or after the photon passes through the shock wave on either 
side of the cut, we use the light-cone perturbation theory |60J to write the SIDIS cross section 
as (see llj for a detailed presentation of a calculation similar to this one, along with the 

appropriate references) 
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where is the well-known light-cone wave function for the 7 * —> qq splitting (normal¬ 

ized as in Q) given by 
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( 86 a) 

( 86 b) 


for transverse (T) and longitudinal (L) polarizations of the virtual photon. Here z = k~ jq~ 
is the light-cone momentum fraction of the photon carried by the tagged quark, cr, cr', A are 
the quark, anti-quark and virtual photon polarizations respectively, i,j are the quark and 
anti-quark colors, and aj = ^(1 — z) Q 2 + m 2 . For simplicity we assume that there is only 
one quark flavor with mass rrif and charge eZf. 

in Eq. (1851) is the fundamental dipole .S’-matrix, with the light-cone Wilson lines 
now, in this low-x regime, stretching to both positive and negative infinities along the x~- 
axis. Note that the interactions of the shock wave with the untagged anti-quark in the first 
diagram of Fig. [ 8 ] cancel on both sides of the cut. 

To extract the unpolarized quark TMD we first impose the Q 2 3>_L 2 condition on Eq. fl85ll . 
with _L denoting any transverse momentum in the problem (this is part of the s> Q 2 3 >_L 2 
condition defining our small -x regime), ft is well-known that the large-Q 2 asymptotics 

re aligned-jet configurations dominated by either 2 —> 1 or z —> 0 
. Since we are interested in producing a quark with k~ & q~ we 


of Eq. (1851) comes from t 
regions of phase space [87 


will only take the z ~ 1 region into account. Note that the transverse photon polarizations 


lead to a leading power of Q 2 . The large-Q 2 limit of Eq. (j85l) with z 
transverse photon polarizations is 
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(87) 


for details of taking the large-Q 2 limit). We have put mq = 0 for simplicity. 
To read off the quark TMD /f we compare Eq. (187|) to the lowest-order SIDIS pro¬ 
cess, modeling the target as a single quark. The dominant contribution to this process is 
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FIG. 9. Lowest-order diagram contributing to the SIDIS cross section on a single quark. 


illustrated in Fig. |9l resulting in the lowest-order SIDIS cross section 

da s T IDIS _ ^ a EM Z 2 


d?kq 




( 88 ) 


Comparing Eq. (1881) to Eq. (1871) and using s ~ Q 2 /x we read off the small-x unpolarized 
quark TMD for a large nucleus to be (see [90] for a similar result) 


where we have also made use of the small-x limit of Eq. (1681) for f ^ uark . 
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(89) 


Having the expression (189j) for the unpolarized quark TMD we can now determine its 
evolution at small-x. It is driven by the small-a; evolution of the dipole S'-inatrix. The 
small -x evolution of S in the leading-logarithmic approximation (LLA) resumming powers 
of a s ln(l /x) is well known and is given by [l 

a, N r f , Q ( x-y ) 2 
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(90) 


with Y = In 1/x ~ ln(s/Q 2 ). In Eq. (l90l) we have separated S xy (Y) = (S) into the Wilson- 
line operator S and the averaging in the nuclear wave function evolved up to rapidity Y (see 
Eq. (fT4]) h Unfortunately Eq. (1901) is not a closed integro-differential equation: the object 
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(, SS ) on its right-hand side is different from (S') on the left-hand side. The equation closes 
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in the large-iV c limit, where it becomes 

a,N r f ,, ( x-y) 2 
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where we again dropped the angle brackets. The initial condition for Eq. (j9Tj) is given by 
(cf. Eq. (HZ])) 
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(92) 


Eqs. (l89j) . m an d (EU) gi ye us expression for and the small -x evolution of the unpolar¬ 
ized quark TMD /j 4 . 



FIG. 10. An example of a diagram contributing a small-x evolution correction to the unpolarized 
gluon TMD. 


The small-x evolution of unpolarized gluon TMD is constructed in a similar manner, with 
some important differences. Instead of DIS with a photon, which only couples to quarks, 
consider “DIS” with a scalar current j = — (l/4)F“ l/ F atiU 


which couples to gluons. A 


gluon dipole diagram contributing to the evolution of the unpolarized gluon TMD is shown 
in Fig. GDI Note that the gluon loop brings in a logarithm of 1/x: hence the diagram is only 
a small-x evolution correction to the gluon analog of the lowest-order graph in Fig. |9] To 
include the small -x evolution into gluon TMDs, one has to start with such a lowest-order 
diagram and first include the Wilson line staple, obtaining the Weizsacker-Williams (WW) 


gluon distribution in the classical approximation 


61 


71 


92 


941 ]. Rewriting the semi-infinite 


adjoint Wilson lines in terms of (derivatives of) the infinite Wilson lines, one can apply 
the JIMWLK evolution obtaining the resulting evolution equation governing the small-a; 
asymptotics of the WW distribution, as was done in !i_7 ! , IS- 

One may worry whether the diagram in Fig. [TO] (and other such corrections) needs to 
be included at large-a:, since the gluon loop may give both logarithms of Q 2 and s. The 
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resolution of this question is in the fact that the energy logarithm coming from the gluon loop 
in Fig. [TO] is ln(s/<5 2 ), which is not large in the s ~ Q 2 regime of Sec. IIII.ll Hence the gluon 
(Fig. m an d quark (Fig. [H]) loop corrections at large-x are suppressed by a power of a s not 
enhanced by a large logarithm, and thus are outside of the precision of our approximation. 

The construction of the small-.x evolution of other TMDs of unpolarized proton or nucleus, 
such as the quark or gluon Boer-Mulders function, can be done similarly to the above, and 
are also governed by the nonlinear evolution of the correlators of fundamental and/or adjoint 
Wilson lines. 

III.2.2. Evolution of Polarized Target TMDs: an Outline 

The evolution of polarized target TMDs at small -x is somewhat different from the un¬ 
polarized case. Here we will only give the general outline of this evolution, with a more 
detailed description left for future work. 

Let us concentrate on the quark TMDs. Again the diagrams of the type shown in Fig. [7] 
need to be considered, except now the quark-gluon vertices are not eikonal, allowing for spin- 
dependence to be transferred from the target nucleus to the produced quark. Note that, as 
a consequence, the contribution of the diagram has a factor of x suppression compared to 
its eikonal contribution to the unpolarized TMDs. These diagrams are supplemented by the 
same-order (both in a s and in x) graphs like that shown in Fig. [TT] 



FIG. If. An example of a leading-order contribution to the polarized-target SIDIS process at 
small-x. 

To include the small-x evolution correction into a generic polarized-target quark TMD 
one has to start by including the effects of GGM/MV multiple rescatterings. The relevant 
diagrams are shown in Fig. [12] There the red bar represents the shock wave again. Note 
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that we are treating the rescattering on a nucleon carrying the spin information about the 
target separately, and draw it explicitly in the graphs. (The spin-dependent scattering may 
also contain gluon exchanges in the t-channel.) Placing such spin-dependent rescattering 
inside the shock wave rectangle, as shown in Fig. [[2] implies that multiple rescatterings may 
happen both before and after the spin-dependent scattering. 



FIG. 12. Diagrams needed for the calculation of the polarized SIDIS cross section and the corre¬ 
sponding TMDs. 

The spin-dependent scattering preserves the z± coordinate of the anti-quark Wilson line. 
Hence all the scatterings of anti-quark and gluon lines at z± cancel between the diagrams 
in the first line of Fig. [12j This is similar to the reason there is no ^-dependence in the 
S x ,y representing the interaction with the target in the first graph of Fig. ( 8 J as shown in 
Eq. (J85J). However, in this case we are interested in the target polarization-dependent TMD, 
and cancellation of the spin-dependent scattering implies that the TMD does not get any 
contribution from the top line in Fig. [T2l We are left only with the two diagrams in the 
second line of Fig. [12] (Note that the photon does not interact with the shock wave it 
crosses: the shock wave to the right of the cut in the first diagram of the second line of 
Fig. [12] only indicates that the 7 * —> qq splitting occurs at positive light-cone time.) 

The same arguments apply to small-x evolution corrections to the diagram in Fig. fTTIfor 
the polarized TMD at hand. In the end one is left with the diagrams where all the evolution 
and interaction with the target happen entirely either to the left or to the right of the cut, as 
illustrated in Fig. [13] For simplicity we only show the linear evolution, realized through the 
quark ladder exchange shown in the left panel of Fig. [T3] and the non-eikonal (non-BFKL) 
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gluon ladder exchange shown in the right panel. The nonlinear evolution corrections are 
represented by interaction with the shock wave. 



FIG. 13. Evolution corrections to the polarizec ■ 
in the right panel is of the non-BFKL type {53 


-target SIDIS process at small-x. The gluon ladder 
. |54|. The red bands represent the shock wave. 


In the left panel of Fig. [I3]the anti-quark in the amplitude emits hard gluons, while itself 
cascading down to lower and lower x until it interacts with the target. After the interaction 
with the target all the emitted gluons recombine back with the anti-quark. The same applies 
for the gluon ladder depicted in the right panel of Fig. [13j Combining each diagram with 
its complex conjugate we write a general polarized-target quark TMD as 

f A 0 , k T ] a, s) = - J - ^ e -*-fe-i) ^2 f ™' A ) 


x 


Ki(y)+Ri:*’(Y) 


(93) 


where A and £ are the quark and the target spin projections correspondingly (taken in either 
longitudinal or transverse basis) and Y ~ In 1/x at small- x. 

Eq. (jU3|) contains two main ingredients. The function f aa i {x — z,y — r, A) has to be sepa¬ 
rately calculated for every polarized TMD (using the 7 * —y qq light-cone wave function) by 
analogy to what was done in the unpolarized target case: this is left for future work. The 
interaction of the x, z (and y, z) dipole with the target, along with its evolution shown explic¬ 
itly in Fig. [13l is included through the Reggeon exchange amplitude R°Z!(Y) {RyZ' *(K)). 
The evolution of (Y) involves mixing between the quark and gluon ladders, as shown in 
the right panel of Fig. [13l It is possible that the evolution of the (non-BFKL) gluon ladder 
depends on the TMD in question. For example, it is conceivable that the evolution of the 
quark transversity could be different from the quark helicity. While further investigation of 
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this evolution for various polarized TMDs is left for future work, let us give an example of 
such evolution by considering only the quark sector. 

By the quark sector we mean considering only the diagram in the left panel of Fig. [[31 
without any mixing with the gluon ladders. The quark ladder evolution is built up out of the 
q —y qG splittings with the quark after the splitting being the softer particle. Corresponding 
small -x evolution kernel is diagonal in quark hclicity. Hence the large- N c nonlinear evolution 
of R^ "' (Y) for a ladder with quarks in the f-channel is independent of the anti-quark helicities 
a, a' on either side of the cut, and in the large- N c limit can be simply read off from the 


evolution of the flavor non-singlet structure function in 

min{y,Y-in[(£-x) 2 /«;2 ]} 


481 ] (see also 


49M551]): 


«y) = r-py) + ^ J^l 


dy rz' + m 


Yi 


(94) 


Here Y) is some initial rapidity for the evolution and r^{Y) is the initial condition. This 
initial condition is given by [48 ] 


x,y 


(Y) = ~ 


s da aa \Y) 


d?b 


exp 


\\x-y\%Ql ln T— 

4 \x-y\ T A\ 


(95) 


with da aa ' (y)/ d 2 b the cross section of lowest-order spin-dependent scattering of the x,y 
dipole on one nucleon in the nucleus (here b = (x + y)/ 2). This cross section has to be 
calculated separately for each TMD. 

To solve Eq. (T94|) one has to first solve Eq. (19T]) to End S'l’|’ 2 00 ’ _00 '(y). Hence the equation 
actually mixes the quark and (BFKL) gluon ladders, albeit at the non-linear level. 

Eq. (19411 is double-logarithmic in energy, that is, it resums powers of a s In 2 s. Hence it con¬ 


tains R , which resums double-logarithms of energy, and S, resumming single logarithms 


48]. 


Indeed, for consistency, the evolution for R has to be augmented by the single-logarithmic 
term: this is left for the future work. Eqs. ([93]) and (f9Tj) complete our outline for the general 
form of the small-.x evolution of the polarized quark TMDs. 

Note again that the full QCD Reggeon evolution (even at the linear level) also includes 
gluon ladders, mixing with the quark ladders, as shown in the right panel of Fig. [131 There 


53 


54|. Hence 


one of the gluons carries the information about the polarization of the target 
Eq. (1941) is incomplete, and has to be augmented by a mixing term with the gluon ladder 
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along with a separate evolution for_the gluon ladder (also mixing with the quark ladder), 


similar to how it was done in |53|, 


54j. While this is outside the scope of this work, here we 


note that the ^-dependence found in [54j scales as 


R aa ~x 


-z s 


(96) 


with the numerical factor z s = 3.45 for 4 quark flavors. The power in Eq. (j96j) is large and 
negative, and can easily become large enough to make the net power of x smaller than — 1 
for the realistic strong coupling of the order of a s = 0.2 — 0.3, resulting in polarized TMDs 
which actually grow with decreasing x fast enough for the integral of the TMDs over the 
low-x region to be (potentially) large. The tantalizing possibility of this growth generating 
a significant contribution to the proton spin coming from low-rc partons will be explored in 
the future work. 


IV. CONCLUSIONS 

In this work, we have presented calculations for the quark TMDs of an unpolarized target 
in the quasi-classical approximation and with leading-logarithmic quantum evolution. At 
face value the calculation applies to a heavy nucleus under the resummation a^A 1 / 3 rs_/ ^(1), 
but its core features should also be valid for a high-energy proton with a large parton 
density generated by quantum evolution. Attempts to model the dense proton as a “nucleus” 
have been successfully applied in the past to phenomenology in inclusive deep inelastic 
scattering, so this approach may be a valuable tool in studying the TMDs of the proton as 
well. Regardless, the formalism presented here makes it possible to perform first-principles 
calculations of TMDs and spin-orbit structure within a controlled resummation of QCD. 

The primary results of our calculation are as follows. The quasi-classical factorization 
formula (1391) expresses the relationship between the quark TMDs of the nucleus, the quark 
TMDs of the nucleons, the distribution of nucleons in the nucleus, and the multiple rescat¬ 
tering on spectator nucleons. The parameterization (l56li expresses the most general form 
of the nucleon Wigner distribution which is consistent with rotational invariance, parity, 
and time reversal. Together, these yield the decomposition of the unpolarized quark dis¬ 
tribution fi (l6Tf and the quark Boer-Mulders distribution h (I78j) in terms of the TMDs 
of the nucleons. The quasi-classical expressions can be taken as initial conditions for the 
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subsequent quantum evolution. In the large-x regime, the leading evolution (l83]l is double- 
logarithmic, resumming powers of a s In 2 Q 2 in agreement with the Sudakov form factor and 
the Collins-Soper-Sterman evolution equation. In the small-x regime, the unpolarized quark 
distribution is given by (1891) . and its evolution is governed by BK-JIMWLK equations (1UU1) 
and (TDTT) . The small-x evolution suggested by the polarized TMDs appears to be more 
complex, involving the familiar BK-JIMWLK evolution as an ingredient for the more intri¬ 
cate Reggeon evolution (1941) (augmented by the mixing with the gluon ladders, as shown in 
Fig. USD. 


Clearly further investigation of the perturbative QCD Reggeon evolution is warranted in 
the future. If the Reggeon evolution leads to polarized TMDs that may grow at small-a; 


faster than 1/x, as Eq. (liMll with the power of x taken from 54] appears to suggest, one 


should try to include saturation effects into this evolution. The inclusion of saturation effects 
may make the integral of the resulting TMD over small -x convergent, and would allow one 
to assess the size of the contribution of such an integral. A large contribution coming from 
the small -x region may help resolve the proton spin puzzle by identifying the phase space 
region containing the missing spin: this possibility is important and has to be explored in 
the future. 


In our formalism, essentially all of the model dependence is encapsulated into the non- 
perturbative Wigner distribution of nucleons inside the nucleus. This distribution, however, 
is highly constrained by symmetry, making it possible to identify distinct channels which 
couple the TMDs of the nucleus to the TMDs of the nucleons. Our approach is also highly 
amenable to modeling and phenomenology. First one chooses an ansatz for the Wigner dis¬ 
tribution such as dnSD; this then determines the form of the quasi-classical multiple scattering 
factor through the integral (l63lh With these ingredients, one can evaluate the intermediate 
integrals over p , b, and r which couple the nuclear TMDs to the nucleonic ones. If desired, 
one can then also choose initial conditions for the nucleonic TMDs, such as the quark target 
model or scalar diquark model of [57], allowing for explicit analytic or numerical calculation 
of the nuclear TMDs. These functional forms can then be evolved with the appropriate 
quantum evolution equations to the kinematics appropriate for comparison with experimen¬ 
tal data. In future work, we would like to extend our computation of TMDs in the dense 
limit to the full set of leading-twist quark and gluon TMDs. This would provide another 
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theoretical benchmark for understanding the range of spin-orbit physics permitted by QCD 
and would be well-suited to a global fit of available TMD data. 


One novel feature of our calculation is that, due to the possibility of spin-orbit coupling in 
the distribution of nucleons, there can be mixing between different TMDs at the level of the 
nucleus and the nucleons. In (Ell) the nucleonic Sivers function mixed into the nuclear quark 
distribution, and in (J78j) the nucleonic transversity and pretzelosity mixed into the nuclear 
Boer-Mulders function. All three of these mixings are due to the same L ■ spin-orbit 
correlation with strength given by Wo am- The fact that the same underlying correlation is 
responsible for multiple mixings is a testable prediction of the theory: in principle, if the 
unpolarized quark distribution ff and Sivers function /of the nucleon are known from 
experiment, then a further measurement of the unpolarized quark distribution ff of the 
nucleus can allow an extraction of the spin-orbit coupling term Woam which is present in 
the nucleus. This would then provide a prediction for the amount of mixing that should occur 
between the nucleonic transversity /if or pretzelosity /i l L f v and the Boer-Mulders function 
h^ A of the nucleus. Such an extraction would require good coverage of the px dependence of 
the nuclear TMDs, and is thus likely to only be accessible at a future electron-ion collider. 


The TMD mixing observed here couples the PT-even and PT-odd sectors, mediated by 
the spin-orbit coupling mechanism and the depth dependence of the multiple scattering. 
This is a very general feature of TMDs in the dense limit which goes beyond the specific 
L ■ Sn coupling that can occur in an unpolarized nucleus. A similar mechanism L ■ Sa was 
already observed for the case of a transversely polarized nucleus, resulting in the mixing of 
the nucleonic unpolarized quark distribution into the nuclear Sivers function 0 . Similar 
features should again occur in the general case of a polarized nucleus, with a relatively small 
number of possible spin-orbit correlations generating a relatively large number of mixings 
between the nuclear and nucleonic TMDs. Indeed, one can imagine performing the same 
sort of analysis for the “generalized transverse-momentum-dependent parton distribution 
functions” (GTMDs) which are analogous operators to ([I]), taken between off-forward matrix 
elements of hadronic states (see, for example, 96]). The GTMDs are the “mother functions” 
from which one can obtain both the TMDs and the generalized parton distributions (GPDs). 
If the same type of mixing occurs at the level of the GTMDs, it raises the interesting 
possibility of a single spin-orbit coupling term being responsible for different mixings in both 
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the TMD and GPD sectors. The methodology presented here, which has at its foundation a 
genuine resummation of QCD, offers a new approach to the calculation of TMDs and related 
quantities, with bountiful applications to both theory and phenomenology. 
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Appendix A: Rotational Invariance on the Light Front 

In Sec. III. 21 we discussed the symmetries of the Wigner distribution of nucleons in the 
nucleus, including the naive expectation that the Wigner distribution should be rotationally 
invariant in the rest frame of the nucleus. Indeed, the Wigner distribution is independent of 
the preferred axis contained in the Wilson lines and virtual photon. But there is a subtlety 
which we will now address: since the Wigner distribution is defined in terms of the light-front 
wave functions, there is a special direction given by the quantization axis with respect to 
which those wave functions are defined. In the standard formulation of light-front perturba- 

so that three-dimensional rotation 
discussed in j^], these light-front 
wave functions are especially well-suited to describing high-energy collision dynamics along 
the 2 -axis, since they are invariant under longitudinal boosts, transverse boosts, and ro¬ 
tations in the transverse plane. But rotations which change the longitudinal direction are 
“dynamical”; that is, they invoke the interaction Hamiltonian and can therefore change the 
particle content of the system. The special role of the 2 -axis also breaks invariance under 
ordinary parity P and time-reversal T, necessitating the application of the modified discrete 
t Note that in this Appendix we will use both c = 1 and c ^ 1 conventions interchangeably: we have taken 
special care to not cause any confusion with this notation. 


tion theory, the system is quantized at fixed x + ~ (ct + z ), 
invariance is broken by the special role of the 2 -axis [j] Ae 
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symmetries “light-cone parity” P\_ and “light-cone time reversal” T± which preserve the 


2 -axis and satisfy P±T± = PT 


98] 


To describe the properties of the light-front wave functions under three-dimensional ro¬ 
tations without invoking the interaction Hamiltonian, it is necessary to use the “covariant 


light-front dynamics” of 


69] in which the quantization axis up is kept arbitrary. If the the¬ 


ory is quantized at fixed uj ■ x, then rotations which transform not only the physical vectors, 
but also the quantization axis uP correspond to “kinematic” transformations which do not 
invoke the interaction Hamiltonian. Consequently, a light-front wave function for spinless 
particles is a Lorentz scalar when the quantization axis a; M is also transformed. In general, 
the spin indices of a light-front wave function rotate according to their group representation, 
but they transform more simply in terms of the covariant spin vector S M 
Defining the lightlike quantization axis as 


^ = VV(1,-") 


(Al) 


for an arbitrary spacelike unit vector n, we see that the analog of the usual “plus” direction 

is 


p + = to ■ p = 



( E p + p-n ), 


(A2) 


where n — z in the usual approach. The Wigner distribution can then be written as 

W xy (p,b-,P;u) = jPy; I <!-«&>•» <JJ((A3) 

where b + = uj ■ b = 0 and the wave functions now explicitly depend on the quantization axis 

Since the “longitudinal direction” is now defined by n, the longitudinal momentum frac¬ 
tion a of a nucleon with momentum p M in a nucleus with momentum is 


a = 


LO ■ p 

U-P 


p+ 

p+’ 


(A4) 


and the vector R M = p 11 — aP dehned so that u ■ R — 0, selects out the momentum 
transverse to n: 


R 2 = (P 0 ) 2 - R 2 


R 2 - {R ■ hf 



(AS) 
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The variables a, R 2 are thus scalars under the generalized transformations which include uP. 


As shown in 


69J, the wave functions depend only on the two invariants: -0 = In 


terms of these quantities, we can rewrite the Wigner distribution (IA3j) as 

W x y(p,b-P-u) = —Lj [ i\SR)'^e-'^ pb e-'^rxia 1 ,K 2 ) ih.(a,R 2 ), (A6) 

J yjaa' 

with 5R p = R p — R' p and 5a = a — a', and project out the polarized and unpolarized 
distributions using (Effij) to obtain 


W(p p , b p , S p - P M ; u p ) = W unp (p p , b p - P M ; uP) - S u W v vol {p p , b p - P **; uP ) 


(AT) 


Thus the Wigner distribution is now manifestly invariant under the generalized Lorentz 
transformations. 

The rotation properties of the wave functions are seen most clearly in the “constituent rest 
frame” (CRF). In light-front perturbation theory, there is no “rest frame” in the conventional 
sense, since one component of the momentum ( p ~) is not conserved. It is thus impossible 
to simultaneously set the net spatial momentum P of the nucleus and JA p t of the nucleons 
to zero. Instead, the momenta satisfy a modified conservation law, 

+ (A8) 

i 

where r is the deviation from the “energy shell.” In the constituent rest frame, one sets the 
net spatial momentum of the constituents to zero, YhiPi — 0 , with the parent particle forced 
to have nonzero momentum, P = Th¬ 
in the constituent rest frame, the variables a,R 2 simplify: 

— R 2 ( = K (p — arn) 2 — [(p ■ h) — ar] 2 = p 2 — {p ■ n) 2 

_UJ-P _ Lo-p _ Ep + p-h c.r.f. E p -\- p • n . . 

a ~ u J -p-^-Pi~E z E i + (E i Pi)-n EiEi ’ 1 J 

where we used (1A8h and u 2 = 0. The wave functions 0(ct, R 2 ) can therefore be expressed 
as 'ijj(p 2 1 p ■ n) which are equivalent invariants but have a simple meaning in the constituent 
rest frame. The covariant formulation thus makes explicit the dependence on the preferred 
direction h; this dependence is a feature of the relativistic wave function and in general 
cannot be avoided. 
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FIG. 14. Simple wave function in which a scalar with mass M and momentum P splits into two 
scalars, one with mass mi and momentum k and the other with mass m2. 


However, there is one limit in which the dependence on the special direction n disappears: 
the non-relativistic limit p 2 <C m 2 or c —> 00. In this limit, the light-front quantization 
condition uj ■ x = ct + (x ■ n) = const, reduces to the equal-time condition ct = const, and 
the dependence on n drops out. Consequently, in the non-relativistic limit the light-front 
wave functions should no longer depend on n: ^(p 2 ,p ■ n) —> ip(p 2 ), restoring the rotational 
invariance of the physical vectors in the rest frame. For example, consider the simple wave 
function for the splitting of massive scalar particles shown in Fig. Q~4l Taking the tri-scalar 
vertex as A, the standard calculation of the wave function gives 


where z 
we have 


^{k,p) = 


1 


A 


2 g+-P + p- 


k — (p — k) 
—A z( 1 — z) 


k + /p A 


(k — zp)? — z( 1 — z)M 2 + (1 — z)m\ + zm% ’ (^-^-0) 

is the momentum fraction of the scalar mi. In the constituent rest frame 


C.R.F. 


i/;(k 2 , k • h) = 


-A z{ 1 — z) 


k 2 — (k ■ n) 2 — 2(1 — z)M 2 + (1 — z)m\ + zm\ 


C.R.F. 

z = 


k 2 + rri 2 + k ■ h 


(All) 


' k 2 + m\ + y k 2 + m 2 

and expanding k 2 -C m 2 , m 2 , M 2 in the nonrelativistic limit (i.e., ( kc) 2 -C (me 2 ) 2 for c —> 00 ), 
we obtain 


-0(P, k ■ n) 


-A 

3 m 2 


1 - 


Ak 2 

3 m 2 


= k 2 ) 


(A12) 
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as desired. Note that in the limit of nonrelativistic constituent motion, the rest frames of the 
parent particle and the constituents coincide, and we no longer need to distinguish between 
the “constituent rest frame” and the ordinary rest frame. 

Indeed, the non-relativistic limit is precisely the limit which is relevant for the orbital 
motion of nucleons in a heavy nucleus. If we expand (IA3D to lowest order in the nucleon 
velocities fd = v/c = \k\/rriN 1, then we have 


E„ « m N (1 + O (I3 2 )) 

asa i(l + IL + 0 (i3 2 )] 

A\ m N V 

P\\ = (P- n) « (a - (1 + O (/?)) 


d 2+ (5p) 

y/p + p ,+ 


d 3 (5p) 

m N 


a +om- 


The Fourier factor exp[— i(p — p') ■ b\ then becomes 


(A13) 


exp[— i(p — p') ■ b] = exp [—ig + _(a — a')P + b ] exp[+z(p — p) ■ 6] 


exp 


—®(P|| ~P\\) [ A7^9+-P + b ) exp[+i(p-tf)-b] (1 + 0 (/?)). (A14) 


Since b + = 0 and b\\ = sj ( b + — b ), we have 

b~ = -V 2 9 + ~ b h 


(A15) 


and in the rest frame 

P +R = \{^(Am N ), (A16) 

such that 

-r —g+-P + b- R = -6,, (A17) 

Attin 

and the Fourier factor becomes 


exp [—i(p — p') ■ b] R = exp [+i(p — p') ■ b\ (1 + O (/?)). (A18) 


Substituting these expressions back into (1A3I) and neglecting the dependence of the light- 
front wave functions on ft for nucleons moving non-relativistically in the heavy nucleus, we 
obtain 


W X y(P } b) R = 


d 3 (p — p') 


2(2tt) 


m N 


e Mp F)-^*(y 2 ) ^ x ,[p 2 ). 


(A19) 
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Once again projecting out the polarized and unpolarized structures using ([26]) . we obtain 


w(p, b, S) R = W unp (p, b) + S • W pol (p, b), 


(A20) 


which recovers the naive expectation of manifest three-dimensional rotation invariance in 
the nuclear rest frame. 

The key is the change of variables between the boost-invariant quantities like a, P + b~ 
and the 3-vectors p, b: 


=>R.F. 

V = 


7* R.F. 

b = 


p, M A \a - — 


b, 


9+- 
' Ma 


P + b~ 


(A21) 


where Ma ~ Am^. Using the dictionary (]A21j) . we can impose rotational invariance in 
the rest frame through p, b to constrain the form of the Wigner distribution. Then we can 
change variables back to the boost-invariant quantities a, P + b~ which describe the high- 
energy collision. 
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